Conservation laws in disordered electron systems: Thermodynamic limit and 

configurational averaging 
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We discuss conservation of probability in noniteracting disordered electron systems. We argue 
that although the norm of the electron wave function is conserved in individual realizations of 
the random potential, we cannot extend this conservation law easily to configurationally averaged 
systems in the thermodynamic limit. A direct generalization of the norm conservation to averaged 
functions is hindered by the existence of localized states breaking translational invariance. Such 
states are elusive to the description with periodic Bloch waves. Mathematically this difficulty is 
manifested through the diffusion pole in the electron-hole irreducible vertex. The pole leads to a 
clash with analyticity of the self-energy, reflecting causality of the theory, when norm conservation 
is enforced by the Ward identity between one- and two-particle averaged Green functions. 

PACS numbers: 72.10.Bg, 72.15.Eb, 72.15.Qm 



I. INTRODUCTION 

Itinerant quantum disordered systems are very dif- 
ficult to describe quantitatively. One has to be aware 
of a few pitfalls that may destroy the standard scheme 
adopted from clean solids. Quantum systems must be 
represented in a Hilbert space. When a static, quenched 
disorder is present, each single configuration of the ran- 
dom potential (force) has its own state space. First as- 
sumption we standardly make is that a suitable represen- 
tation space exists comprising states of all configurations 
of the random potential. The disordered system then for- 
mally behaves as a clean one with renormalized products 
of operators via vertex corrections. 

Operators are generally difficult to treat unless we are 
able to use their spectral representations. In the thermo- 
dynamic limit, i. e. in infinite volume, the eigenstates of 
itinerant Hamiltonians may either be localized, quadratic 
integrable, or extended with a plane-wave character. The 
former eigenfunctions belong to the point spectrum with 
separate eigenvalues, while the latter ones to a continuous 
band of eigenenergies. These two types of eigenfunctions 
are well separable for one inhomogeneous configuration. 
Configurational averaging, however, can turn the point 
spectrum effectively continuous and it becomes extremely 
difficult to distinguish the character of spectral energies, 
that is the extension of the eigenfunctions, solely from 
averaged quantities. Determination of the extension of 
typical eigenstates of a random itinerant Hamiltonian is 
the core of the so-called Anderson localization problem. 1 

Since exact solutions of the Anderson localization 
problem are not available one has to resort approximate 
treatments. There are two complementary tools, numer- 
ical and analytic, to tackle this problem. Neither of them 
is, however, able to answer all questions about the spec- 
trum of random itinerant Hamiltonians. While the for- 
mer technique can deal only with finite lattices and many 
configurations of the random potential^ the latter one 
deals mostly with the thermodynamic limit and config- 



urationally averaged quantities^ In particular the latter 
approach demands a special caution in building a con- 
sistent picture of the metal-insulator transition at which 
configurations with extended states become irrelevant. 

In order not to miss the physics of the problem one 
has to minimize spurious effects, introduced in approxi- 
mate treatments, by taking into account as many exactly 
valid relations as possible. Conservation laws are typical 
properties to be maximally satisfied also in approximate 
theories. They impose the so-called "no-go" theorems 
that may essentially restrict the state space of the sys- 
tem and hint to available "order parameter" . The aim 
of this paper is to elucidate the role and applicability of 
conservation of probability, or conservation of the norm 
of the wave function, in the Anderson localization prob- 
lem. Although it may seem evident that the norm of the 
wave function must be unconditionally conserved, it is 
the random character of the Anderson model that makes 
the proper utilization of probability conservation a rather 
subtle issue deserving a deeper and more thorough inves- 
tigation. 

We use the simplest possible model on which we can 
demonstrate subtleties of the application of conservation 
laws in random systems. It is the Anderson model of 
noninteracting lattice electrons subject to a random po- 
tential, the Hamiltonian of which can be written as 



Had = *« C < C J + H ViC i Ci ■ 



(1) 



We assume that the kinetic energy (nearest-neighbor 
hopping term) tij — t\i-j\ is homogeneous and only the 
local potential Vi is random. The values of the random 
potential Vi are site-independent and obey a disorder dis- 
tribution p(V). A function depending on the random 
potential Vi is then averaged via 



<*(V5)>o, 



dVp(V)X(V) 



(2) 
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Anderson showed in his seminal paper, Ref. [l, that 
the eigenstates of single configurations can either be ex- 
tended or localized. In the former case we have a metallic 
conductivity (diffusive regime). The latter case, called 
Anderson localization, leads to an insulator and we have 
a diffusionless transport regime. A consistent theory of 
the disorder induced vanishing of diffusion is still miss- 
ing. One of the reasons for this shortage is incomplete 
understanding of probability conservation in this simple 
model. We discuss in the following sections what holds 
for sure and what applications of conservation laws de- 
mand special additional assumptions. 



II. NORM CONSERVATION AND WARD 
IDENTITY: ONE CONFIGURATION 



We start with a finite system and consider one fixed 
configuration of the random potential. It is a very 
straightforward situation where no problems with con- 
servation of the norm of the wave function are to be ex- 
pected. Since the Anderson model, Eq. (JTJ, consists of 
noninteracting electrons, we can formally find an exact 
solution via an evolution operator. We define 



U(t) 



-iHt 



(3) 



The evolution operator U(t) represents a formal solution 
unless we are able to diagonalize the Anderson Hamil- 
tonian. It is, however, possible only numerically and for 
rather small lattices. We can nevertheless use this formal 
solution to see how conservation of the norm of the wave 
function can be represented. We evaluate the following 
integral for arbitrary complex energies z with 3z > 



dt 



(m\m) 



dt e izt 



1> 



dt A i>*j 



* izt , 

e {j 



U(t) 



U{t)U\t) 

7 y £ ^ 



where \i) denotes the localized orbital on the site and 
<*|j> = %• 

Identity J3J is a rather evident consequence of conser- 
vation of the norm of the wave function. It is less evident 
and more interesting when we rewrite this identity with 
resolvent operators and Green functions. We introduce 
the resolvent as a Fourier transform of the evolution op- 
erator 



G(z) = 



dt 9(t) e lzt U{t) = 



zl-t-V 



(5) 



Using the local orbitals (Wannier states) we define one- 



and two-particle Green functions 



G ij {z) = (i\d{z)\j) 



G{z 2 ) 



(6a) 
(6b) 



The two-particle Green function is trivially only a prod- 
uct of one-electron Green functions, since we do not have 
interacting electrons and we stay within one fixed config- 
uration of the random potential. 

It is easy to insert the resolvent operators (finite ma- 
trices) into Eq. (@J to obtain 

E G Snj( Zl ' 22 ) = ^G in {zi)G n j{z 2 ) 

n n 

l 



z 2 - Zl 



[Gijiz!) - G l3 (z 2 )} . (7) 



The result is kind of Ward identity being a direct con- 
sequence of probability conservation for quantum nonin- 
teracting systems. It holds for any configuration of the 
random potential. 

In finite systems all eigenenergies of the Anderson 
Hamiltonian are discrete and nondegenerate. The sit- 
uation changes if we perform the thermodynamic limit, 
where the volume of the lattice increases to infinity, i. e., 
we have infinite-many lattice sites. To perform the ther- 
modynamic limit we choose the values of the random 
potential on the newly added sites according to the dis- 
order distribution p from Eq. (0). To extend the validity 
of the Ward identity (JJJ) to infinite systems we have to 
demand completeness of the set of Wannier states. That 
is, we require 



1 



E (tp \n) (n| ip) 



1 



(8) 



The latter equality says that we close finite sums over the 
Wannier states in such a way that the norm can be fixed. 
The Wannier states form an orthonormal basis in the 
Hilbert space of random configurations, i. e. — Sij. 
Condition JHJ guarantees validity of the Ward identity J2J 
strictly only for proper states from the Hilbert space (lo- 
calized states). In this case the infinite sum on the l.h.s. 
of Eq. (0 converges strongly in the Hilbert space spanned 
over the Wannier states. We, however, know that the 
thermodynamic limit generally generates extended states 
belonging to the continuous spectrum of the Anderson 
Hamiltonian. These extended states are not proper vec- 
tors in the Hilbert space. To prove the Ward identity also 
for extended states of a single configuration we have to 
weaken the convergence criterion for the sum over the in- 
termediate states. It is in principle possible to formalize 
such a procedure. 



III. CONFIGURATIONALLY AVERAGED 
QUANTITIES 

Analysis of individual configurations of the random 
potential is manageable only by numerical means, since 
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random configurations are not generally translationally 
invariant. Hence, an analytic treatment is not viable due 
to the lack of apparent symmetry that would help reduce 
the number of degrees of freedom. 

The configurational averaging provides a means for 
reducing the complexity of the description of random 
systems, as it restores translational invariance. Trans- 
lationally invariant lattice electrons can be formally de- 
scribed by standard field-theoretic methods used for pure 
(interacting) systems. If the ergodic (self-averaging) hy- 
pothesis holds, i. e., if extensive thermodynamic quan- 
tities evaluated for one typical configuration equal their 
averaged values, averaging does not change the physics 
of the problem. It does not, however, mean that all re- 
lations and identities valid for single configurations can 
automatically be transferred to averaged quantities. 

To generalize the Ward identity Q to averaged quan- 
tities we have to assume that the sum over the intermedi- 
ate states commutes with the configurational averaging, 
that is 

^E*«-E^E*»(M), o) 

C {C} n n c {C} 

where \1/„({C}) are configurationally dependent func- 
tions of interest, the matrix elements of the two-particle 
Green function in this case. This equality may seem evi- 
dent and indeed holds for every finite- volume samples as 
all involved sums run over finite number of elements. In 
the thermodynamic limit both sums, over the basis states 
and over the random configurations Af^ 1 X){c}> are 
infinite and they are interchangeable only if states of all 
relevant random configurations fall into a single Hilbcrt 
space. That is, if the basis vectors \n) are configura- 
tionally independent. We will show later on that this 
need not be always the case. 

After averaging the system does not behave as a sim- 
ple Fermi gas, since the product of averages does not 
equal the average of the product. Hence, the two-electron 
resolvent, Eq. Q6bfl. will be defined via a Bethe-Salpeter 
equation, 

(G^ kl (zi, z 2 )) = (G«(*i)> „ (G M {z 2 )) av 
+ £ (Gn'(zi)) av (Gi'i(z 2 )) av 

i' j' k'V 

x Aiwfazi) ( G ?lkk'( z ^ z ^) av ■ ( 10 ) 

We introduced a two-particle irreducible vertex A to take 
into account correlations between the one-electron Green 
functions from the product in Eq. (|6b() due to the averag- 
ing. Analogously, the averaged one-electron resolvent JSJ 
is determined from a self-energy X via the Dyson equa- 
tion, (C?(z)) = [zl — t — S] . Both, the self-energy 
£ and the irreducible vertex A, display translational in- 
variance. 

Assuming validity of Eq. © and using the Bethe- 
Salpeter equation (fTUf> we transform the Ward identity 



to an identity between the self-energy and the two- 
particle irreducible vertex^ 

Eu(zi) - Y>u{z 2 ) 

= ^2Mj',k'i(zi,z 2 ) [(Gj>k'(zi)} av ~ {Gj>k'(z2)} av \ ■ 
j'k' 

(11) 

It is a special case of the so-called Vollhardt-Wolfle Ward 
identity proved from a diagrammatic perturbation expan- 
sion for the two-particle vertex A£ 

To exploit fully the simplifications brought in by the 
configurational averaging we have to switch to a basis set 
consisting of Bloch waves. It is a natural basis for trans- 
lationally invariant lattice systems. The Bloch waves can 
be introduced through the Wannier states as follows 

(k|n) = e -* k - R - , (12) 

where k is quasiparticle momentum labelling the Bloch 
waves. As far as we stay in finite volumes, the complete- 
ness of the basis and the norm of any state vector read 

l^|k)(k|=T, IMI 2 ^ $>;|k)(k|^} (is) 

k k 

with N being the number of sites in the given volume. 
The only difference compared to relations JSJ is the nor- 
malization factor of the basis vectors, (k|k') = A(5k,k'- It 
says that the Bloch waves are normalized to the volume, 
i. e., their norm increases linearly with the volume. This 
normalization is chosen so that the thermodynamic limit 
is well defined. There the completeness relation comes 
over to a spectral integral 

-J-V|k)(k| > -i-, f d d k\k) (k| =T (14) 

N ^ N-,oo (27r) d J 1 x v ; 

with generalized orthonormality (k|k') = <5(k— k'), where 
5 is the Dirac delta function. 

In the Bloch-wave representation the one-particle 
Green function is diagonal 

G(k ' Z)= z-e(k)-S(k,z) 

4E r!MR,_Rj> ^WL ( 15 ) 

and the two-particle Green function depends on only 
three momenta as the total momentum is conserved. The 
notation for the momentum-dependent averaged two- 
particle Green function we use throughout the rest of 
the paper reads 

ijkl 

x^^e^jG^fz,^)) . (16) 
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The Ward identity Q for the averaged Green functions 
in momentum representation is 

E 4% (^1^2; 0) = - J- [G(k,zi) - G(k, z 2 )] 

k' 

(17) 

where Az = z\ — z 2 - 

Notice that in finite volumes one cannot straightfor- 
wardly distinguish between diffusive electron states, de- 
scribing a particle "smeared" throughout the whole sam- 
ple, and localized states with exponential tails confined 
to a finite subvolume. Such a distinction is meaning- 
ful only in the thermodynamic limit. There the diffusive 
states have a character of Bloch waves, i.e., they describe 
particle fluxes. On the other hand the localized states 
contain just a single particle localized in a finite volume. 
Intuitively, technical difficulties can arise if a physical 
process involves both types of eigenstates, localized and 
extended. It is the case of the Anderson localization tran- 
sition where extended states go over to localized ones. 

IV. DIFFUSION POLE AND CAUSALITY 



function with uncorrelated part subtracted, 

(z + , z_ ; q) = G(k, z+)G(k + q, z_) [S(k - k') 
+r kk ,(z+, z_; q)G(k', z+)G(k' + q, z_)] . (20) 

In analogy to the averaged two-particle resolvent, the 
two-particle vertex obeys a Bethe-Salpeter equation. 
Once we go beyond local approximations, an example 
of which is the mean-field coherent-potential approxima- 
tion (CPA), there are three Bethe-Salpeter equations for 
the vertex T we can construct. The three Bethe-Salpeter 
equations are related to three different ways how to de- 
fine a two-particle irreducibility. Equivalently they cor- 
respond to three topologically distinct two-particle pro- 
cesses, the electron-hole and the electron-electron (hole- 
hole) scatterings and the one-electron self-corrections^ 

Since only the off-diagonal one-electron propagators 
G(k, z) = G(k, z) — N^ 1 ^ k G(k, z) are relevant for dis- 
tinguishing various representations of the full vertex we 
use Bethe-Salpeter equations with just these propaga- 
tors. The most important Bethe-Salpeter equation is 
that from the electron-hole scattering channel that reads 



Since the averaged two-particle resolvent is not a sim- 
ple product as for single configurations, Eq. (J5J, but 
rather fulfills the Bethe-Salpeter equation . it is not 
so apparent how the Ward identity i|17|) is related with 
particle conservation. It can be demonstrated at best via 
diffusion. In the weak disorder limit the motion of elec- 
trons in a random medium is expected to be diffusive, 
i.e., the Fourier components of non-equilibrium particle- 
density variations die out according to the exponential 
law 



5n(t, q) = Sn(0, q) e 



-Dq 2 t 



t>0. 



(18) 



In order to reproduce such a diffusive behavior 
the electron- hole correlation function <5^(q, u) — 

N ~ 2 Ekk' G kk' ( e f + w + i0+, E F - i0+ ; q) has to display 
the so-called diffusion polep 



2nnp 



-iio + Dq 2 



(19) 



Here D stands for the diffusion constant and np is the 
electron density of states at the Fermi energy Ep. The 
diffusion pole exists only if the Ward identity l|17|) is 
fulfilled by the averaged one- and two-particle Green 
functions^ 

Unfortunately, the correlation function $ does not 
obey an equation of motion, so it is difficult to use this 
realization of the diffusion pole to control conservation 
laws in approximate theories. It can, however, be shown 
that the same pole must appear in the two-particle irre- 
ducible function A. To demonstrate that we analyze the 
two-particle functions in more details. It is convenient to 
introduce a vertex T that is just the two-particle Green 



rkk'(q)=A k ' k ,(q) + l^A k ' k „(q) 



N 

k" 

X G+(k")G_(k" +q)r k „ k ,(q) . (21a) 

We added the superscript eh to the irreducible vertex A 
in order to distinguish it from the other two irreducible 
channels. Notice that the full vertex T must remain chan- 
nel independent. For the sake of simplicity we suppressed 
energy variables. The subscript ± at the one-electron 
propagators refers to the first and second energy variable 
in the corresponding two-particle vertex, respectively. 

Besides the electron- hole channel, Eq. I|21a|) . it is 
for our purposes sufficient to use only one out of the 
two remaining Bethe-Salpeter equations — the electron- 
electron (hole-hole) channel, 

T kk , (q) = A-, (q) + 1 £ A&,, (q + k' k") 
k" 

xG + (k")G_(Q-k")r w (q + k-k") , (21b) 

where Q = k + k' + q. 

The two-particle irreducible vertices play the role of 
the two-particle self-energy and are to be determined 
from e. g. a diagrammatic expansion. The best way 
to access them is to use the so-called parquet approach,- 
In it we introduce a completely irreducible (irreducible 
in all distinct channels) two-particle vertex X and use 
topological nonequivalence of different scattering chan- 
nels. When we distinguish only the electron-hole and 
electron-electron channels we can write the parquet equa- 
tion as^ 



IW(q) = Af k ,(q) + A k e k ,(q) - Z kk ,(q) . 



(22) 



■5 



The minus sign at the completely irreducible vertex X 
compensates for the same contributions in the electron- 
hole and the electron-electron irreducible vertices. No- 
tice that all local two-particle contributions belong to 
the completely irreducible vertex X. 

If the physical system under investigation is invari- 
ant with respect to time reversal (no magnetic field), the 
two-particle vertices obey the electron-hole symmetry ex- 
pressed in a relation Tkk' (q) = Tkk' (— q — k — k') for the 
full vertex and in Agj c ,(q) = A^, (— q — k — k') for the ir- 
reducible ones. As a consequence of the parquet equation 
<|22|) the diffusion pole from the full vertex T must appear 
at least in one of the irreducible functions A. However, 
due to the electron-hole symmetry, both vertices A eh and 
A ee display the same analytic behavior and hence must 
contain the same singularity, the diffusion pole. As dis- 
cussed above, the existence of the diffusion pole in the 
two-particle irreducible vertices is a consequence of con- 
servation laws applied to averaged systems via the Ward 
identity. Since the vertex functions obey equations of 
motion, the existence of the diffusion pole in them can 
be used to check conserving character of approximate 
treatments. 

Having a pole in the electron-hole irreducible vertex 
A eh , however, may have unexpected implications. The 
electron-hole irreducible vertex is bound with the self- 
energy via the Ward identity. The singular behavior of 
the electron-hole irreducible vertex could then be trans- 
ferred onto the self-energy. This can be made transparent 
by evaluating the function 



AW(uj) 



1 

N 



E - lu + -Y,(k, E + uj + i0 + )} 

(23) 

in such a way that the self-energies are expressed in terms 
of the two-particle vertex A eh according to the Fourier 
transform of Eq. (jlip. Doing so we finally find out£ that, 
besides an analytic part, the function AW contains also 
a singular term of the form 



AW s d m9 {io)^K\n 2 F 









Dk% 


u 


OJ 




Dk% 



d/2 
d/2 



In 



Dki 



for d ^ 41, 
for d = 41, 



(24) 



where K is a dimensionless constant, kp stands for the 
Fermi momentum and d denotes spatial dimension. The 
non-analyticity of the self-energy, and therefore of the 
one-electron Green function, is unacceptable as it con- 
tradicts another very fundamental property — causal- 
ity. Causality is a consequence of self-adjointness of the 
Hamiltonian, that is of reality of eigenenergies. It can- 
not be given up in any physically meaningful treatment. 
We are thus led to a surprising conclusion that the Ward 
identity cannot survive to the thermodynamic limit 
of configurationally averaged Green functions. At least 
for real-energy differences denoted ui in Eq. (|23|) and (12411 . 
This indicates that there must be a flaw in the seemingly 
straightforward extension of the Ward identity from fi- 
nite to infinite configurationally averaged systems. 



APPROXIMATE ANALYTIC TREATMENTS 



It is impossible to resolve rigorously the incompatibil- 
ity of probability conservation and causality in infinite, 
configurationally averaged systems, since exact solutions 
are not provided. We hence have to resort to approx- 
imate treatments to trace down the problem and clues 
for a possible resolution. 

It is natural to start with a mean-field solution for 
disordered electron systems, being the exact solution of 
the problem in infinite spatial dimensions. This limit 
amounts to the coherent-potential approximation. We 
have explicit expressions for all quantities in this ap- 
proximation. The self-energy is determined from a Soven 
equation that can be cast into 



G(z) = ([G-\z) + nz)-V l Y 1 



(25a) 



where G(z) = N 1 J2k C(k, z). The two-particle irre- 
ducible vertex then is 



A(z+,z-) = 



f 



G(z+)G(z.) 



1 - 



f 



1 



1 + (£(z+) - Vi) G(z+) 1 + (E(z_) - Vi) G(zJ) 



S(z+)-S(z_) 
G{z+)-G{z-) 



(25b) 



r 



The self-energy from the Soven equation (|25a|l can be 
proven analytic and the last equality in Eq. (|25b|l ex- 
presses the Ward identity. The conservation law is ful- 



filled and apparently does not contradict causality. 

The CPA is indeed a consistent theory, but only if 
we keep strictly to local quantities. To remain consis- 
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tent at the two-particle level we should hence consider 
only the local vertex being 7(2+, zJ) — \(z+,Z—)/(l — 
A(z + , z_)G(z + )G(z_)). This local vertex, however, does 
not contain the diffusion pole. Velicky derived in Ref. [U 
the full two-particle vertex in the CPA in a form 



r kk '(q) = 



A 



l-Ax(q) 
where x is a two-electron bubble 

^) = ^E G +( k ) G -( k+c i) 



(26) 



(27) 



The CPA nonlocal two-particle vertex, Eq. I|26|) , contains 
the diffusion pole precisely as in Eq. I|19f) . see Ref. for 
explicit calculation. Even at this level the CPA seems 
to be consistent. The CPA does not treat, however, 
the electron-hole and the electron-electron vertex equiv- 
alently. We have A eh — A, whereby there is no expres- 
sion for the electron-electron irreducible vertex A ee . The 
clcctron-holc symmetry for two-particle vertices is there- 
fore broken in the mean-field solution. As a consequence 
of the violated electron-hole symmetry at the two-particle 
level we have no sign of localized states in the CPA. 

To include localized states and to approach the An- 
derson metal-insulator transition we have to go beyond 
local approximations. Already the notion of weak lo- 
calization needs a sum of infinite-many "crossed" two- 
particle diagrams (ladders from the electron-electron 
channel) leading to a singular electron-hole irreducible 
vertex. The importance of the electron-hole symme- 
try at the two-particle level and the diffusion pole in 
the electron-hole irreducible vertex was stressed and ex- 
ploited by Vollhardt and Wolfle in their theory of Ander- 
son localization^ 

Theory of Vollhardt and Wolfle is a self-consistent ap- 
proximation for the electron-hole and electron-electron 
irreducible vertices. Due to the complete electron-hole 
symmetry we obtain a nonlinear equation for a generic 
vertex. This vertex displays the diffusion pole of Eq. (|19(l . 
Further on, the Ward identity l(TT)l is assumed to hold. 
In this situation the two-particle irreducible vertex in the 
hydrodynamic limit (q — > 0) contains only a single free 
adjustable parameter, the diffusion constant D. Ander- 
son localization is signalled in this approach by vanishing 
of the diffusion constant D. 

A diagrammatic expansion for two-particle irreducible 
functions is formally used to select an appropriate ap- 
proximation in the Vollhardt- Wolfle theory. This the- 
ory, however, does not provide a similarly closed expres- 
sion for the one-electron self-energy. Since the spectral 
behavior of the self-energy is not decisive for the exis- 
tence of the Anderson localization transition, the self- 
energy in the one-electron propagators is taken from the 
weak-scattering limit, i. e. from the Born approximation. 
The Ward identity is used only to recover the diffusion 
pole, but not to determine the self-energy consistently 
so that it be compatible with the chosen approximate 



vertex function. It is then evident that we cannot check 
compatibility of the Ward identity with causality within 
the Vollhardt- Wolfle theory. Due to Eq. (|2IJl we cannot 
expect that a causal self-energy could be found to the 
approximate vertex resulting from the Vollhardt- Wolfle 
approach. 

We can use another route to trace conservation laws 
in configurationally averaged infinite systems. Since the 
CPA obeys all restrictive conditions except for the two- 
particle electron-hole symmetry, we can try to correct 
this approximation just by implementing this missed fea- 
ture. This can actually be achieved within the par- 
quet approach providing a framework to sum systemati- 
cally nonlocal vertex corrections to the CPA irreducible 
vertex^ We showed recently how to solve the parquet 
equations in electron-hole symmetric theories asymptot- 
ically exactly in the limit of high spatial dimensions. 9 
The solution for the full two-particle vertex resembles 
that from the CPA, but the electron-hole symmetry is 
manifestly present. We obtain 9 



r kk '(q) = 7 



Ar 



A x(k + k' + q) 



Aox(q) ^ 

1 - A oX (q) ^ 1 - A x(k + k' + q) 



(28) 



where we used the local CPA vertex 7 and denoted 
A = A /(1 + A G + G_) andx(q) =x(q)~G + G_. The 
only local parameter to be determined from the parquet 
equation is An. It fulfills an equation 



A = 7 + A 



1 x ^ Anx(q) 
N^l- Aox(q) 



(29) 



It is clear that the CPA vertex, Eq. I|26|) . is reproduced 
from Eq. (|28|l if we put An = 7 in it and if the last term 
on its right-hand side is neglected. It is just the neglected 
term that restores the two-particle electron-hole symme- 
try. On the other hand, it is the same term that makes it 
difficult to comply with the Ward identity. The vertex T 
from Eq. (|28|l does not obey a Bethe-Salpeter equation. 
It is a sum of two solutions of the Bethe-Salpeter equa- 
tions in the electron-hole and the electron-electron chan- 
nels. It is hence unclear how to choose the self-energy so 
that it would fit the Ward identity. 

To restore the diffusive behavior of the full vertex we 
have to use the Ward identity for a specific combination 
of energies and determine the self-energy from the irre- 
ducible vertex An. We use 9 - 

3£(£ + irf) = A {E + ir], E - ir))SG(E + irj) . (30) 

to resolve the imaginary part of the self-energy. The 
real part of the self-energy is then calculated from the 
Kramers-Kronig relation 



dE' + irj) 



E' - E 



(31) 



In this way we construct a consistent approximation 
with the two-particle electron-hole symmetry and with 
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a causal self-energy that satisfies the Ward identity in 
a maximally possible way. Since the Ward identity was 
used in Eq. (|3U|) only for imaginary energy difference in 
the vertex function, this solution is not in conflict with 
Eq. (|24|l where the energy difference is real. 

The Ward identity does not hold in this solution for 
real energy differences. As a consequence the diffusion 
pole from Eq. (|f 9fl modifies to 



2irnp 



-iA{w)u + D(u)q 2 



(32) 



The weight of the pole is no longer one, A(Q) > f. 
The quasistatic and hydrodynamic limit of the vertex 
function is then governed by two parameters, the diffu- 
sion constant D(0) and the weight of the diffusion pole 
(wave- function normalization) A(O)" 1 . The Anderson 
localization transition is signalled here by vanishing of 
the weight of the diffusion pole j4(0) -1 — > accom- 
panied by vanishing of the effective diffusion constant 
D' = D(0)/A(0) -> 0. 



VI. PROBLEMS WITH A QUANTITATIVE 
DESCRIPTION OF LOCALIZED STATES IN 
TRANSLATIONALLY INVARIANT SYSTEMS 

Our analysis disclosed a surprising fact that it is im- 
possible to describe consistently disordered electron sys- 
tems with averaged Green functions so that conservation 
of particle number would hold in the state space of Bloch 
waves, the eigenstates of the pure solids. Since the prob- 
ability density is not conserved in the description with 
averaged functions, we have to ask whereto can particles 
vanish? 

To answer this question and to understand deviations 
from the Ward identity in averaged infinite systems, we 
recall that there are two very different types of states 
in the thermodynamic limit — extended and localized. 
The former are normalized to volume. Their normal- 
ization factor increases linearly with increasing volume. 
The latter are normalized to unity. Their normalization 
constant is effectively independent of the volume of large 
samples. This difference may seem to be a formal issue, 
but in fact it poses a principal obstacle for describing 
both types of states within one theoretical framework. 

Each configuration of the random potential in a disor- 
dered system has its own representation space. It is the 
Hilbert space spanned over the eigenstates of the Hamil- 
tonian with given values of the random potential. We 
can generally characterize the representation space via a 
decomposition of unity operator 



1 



(33) 



where we distinguished extended eigenstates |x) and lo- 
calized ones \u). The total number of eigenstates is 



Next + Ni oc — N. To determine the eigevalues of the con- 
figurationally dependent Hamiltonian we have to spec- 
ify boundary conditions for the eigenvalue problem. We 
standardly use periodic boundary conditions for the ex- 
tended states and asymptotic vanishing of the localized 
states. It is evident that we cannot describe localized 
states with periodic boundary conditions, unless they are 
periodically repeated as are e. g. Wannier states. Con- 
figurationally dependent localized states are confined to 
a specific finite part of the space and do not replicate 
themselves periodically, since one configuration is typi- 
cally not translationally invariant. 

We cannot effectively work with individual configu- 
rations, since we know neither the extended nor the lo- 
calized eigenstates. Moreover, the ratio Ni oc /N ext varies 
from configuration to configuration. In the thermody- 
namic limit we are able to deal only with extended 
or periodically distributed localized states. We stan- 
dardly use the Hilbert space spanned over the Bloch 
waves and assume that all relevant configurations of 
the random Hamiltonian can be accommodated in this 
space. There are no problems with delocalized eigen- 
states, since we have (x|k)/(k|k) = O(l). However, the 
localized, translationally noninvariant states fall out from 
this Hilbert space, since (z/|k)/(k|k) = OlN^ 1 ). Hence, 
if we have macroscopically many configurations with lo- 
calized eigenstates, we observe deviations from the Ward 
identity for the averaged Green functions and probability 
will no longer be conserved in the space of asymptotically 
free states. Mathematically it means that the number of 
extended eigenstates N ext is configurationally dependent 
and we cannot interchange the configurational averag- 
ing with the summation over the intermediate states in 
Eq. ©. 

Evasion of localized states from the space of extended 
waves is a severe drawback of the description, since we 
loose the normalization of the density of states. To repair 
this we replace the missing localized states with extended 
ones so that we keep the number of available states always 
equal the number of the lattice sites N . Thereby we 
believe that the differences between the exact unknown 
localized and the supplemented known Bloch waves can be 
represented by corrections via the self-energy and other 
irreducible vertex functions. We recover the sum rule 
for the density of states, but the Ward identity at best 
only in the quasistatic limit lu — > 0. This is the price we 
pay for our inability to describe the localized states in 
translationally invariant systems exactly. 



VII. CONCLUSIONS 

We discussed in this paper applicability of conserva- 
tion laws in disordered itinerant noninteracting systems. 
We showed that formal manipulations leading to conser- 
vation of the norm of the wave function are justifiable 
only for individual configurations of the random poten- 
tial in finite volumes. Extensions of the conservation 
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laws via Ward identities to infinite systems demand a 
caution when dealing with infinite sums. In the thermo- 
dynamic limit we have to distinguish proper states from 
the Hilbert space with a finite norm (localized states) and 
generalized states with the norm (or the normalization 
constant) proportional to the volume (extended states). 
If we stay within individual configurations of the random 
potential we can well separate these two types of states 
and rely on the derivation of the Ward identities. 

Problems arise, however, when we start to average 
over configurations to restore translational invariance. A 
generalization of the Ward identity to the averaged Green 
functions faces principal obstacles, when both types of 
states, localized and extended, are present. An exten- 
sion of the Ward identity to the averaged functions can 
work when eigenstates of all relevant configurations are 
from the same Hilbert space. Based on the observed be- 
havior of averaged systems we conclude that, whenever 
localized states become macroscopically relevant, we are 



unable to find the proper Hilbert space to represent the 
localized eigenstates exactly. We have at our disposal 
practically only the Hilbert space spanned over the Bloch 
waves. The localized states are in the thermodynamic 
limit orthogonal to Bloch waves, (f|k)/ (k|k) = 0(N^ 1 ), 
and hence elusive to a quantitative description within the 
available Hilbert space. Since the space of Bloch waves 
is incomplete, probability is no longer conserved, when 
macroscopic portion of configurations contain localized 
states. Our standard assumption that configurational 
averaging does not influence the representation space 
and can be represented as averaging over rotations in 
the given Hilbert space is hence incorrect when quantum 
coherence between spatially distinct scattering events is 
taken into account. 
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